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We investigate the dynamics of a magnetic system consisting of two magnetic moments coupled 
by either exchange, dipole-dipole, or Dzyalozhinski-Moriya interaction. We compare the switch- 
ing mechanisms and switching rates as induced by the three couplings. For each coupling and 
each configuration of the two anisotropy axes, we describe the switching modes and, using the 
kinetic theory of Langer, we provide (semi-)analytical expressions for the switching rate. We then 
compare the three interactions with regard to their efficiency in the reversal of the net magnetic 
moment of the dimer. We also investigate how the energy barriers vary with the coupling. For the 
dipole-dipole interaction we find that the energy barrier may either increase or decrease with the 
coupling depending on whether the latter is weak or strong. Finally, upon comparing the various 
switching rates, we find that the dipole-dipole coupling leads to the slowest magnetic dimer, as 
far as the switching of its net magnetic moment is concerned. 



I. INTRODUCTION 



Multilayered magnetic systems as permanent mag- 
nets with high performances encounter a renewed in- 
terest due to their potential use as magnetic record- 
ing media with high thermal stability and reduced 
switching fields^ ^. Another candidate for informa- 
tion storing media, also of constantly growing in- 
terest, is provided by magnetic nanoparticles. For 
such applications, it is very important to understand 
how the dynamics of these systems and their reversal 
mechanisms are altered by a change in their physical 
parameters such as the underlying material, thick- 
ness, and stacking conditions. There are also rather 
involved issues related with surface and interface ef- 
fects even if one assumes, to some extent, that the 
problems of crystalline diffusion are properly dealt 
with during fabrication. Apart from the issues re- 
lated with the intrinsic properties of the constituting 
elements (layers or particles), an issue of paramount 
importance is that of (inter-layer or inter-particle) in- 
teractions as these affect energy barriers which they 
have to circumvent during their switching process. 
Now, because the latter play a central role in mag- 
netic recording technology, the problem of interac- 
tions must be addressed in a broader way. Indeed, 
several types of interactions may occur in particle en- 
sembles or in multi-layered systems. However, it is a 
difficult task, if at all possible, to tell in detail what 
interactions are involved in the dynamics of the sys- 
tems to be studied. In general, one resorts to the 
"molecular field" approach and considers only the ef- 
fective interaction. On the other hand, it is clear that 
the dynamics strongly depends on the type of interac- 
tion that is operating within the system. As such, it 
would be useful to relate the observed dynamics to a 
given type of interaction, even as an effective one. In 
spin systems, there are short-range as well as long- 
range, isotropic or anisotropic, interactions. These 
are mostly exchange interactions (EI), dipole-dipole 



interactions (DDI), or Dzyalozhinski-Moriya interac- 
tions (DMI). 

On the atomic level, for instance, these three inter- 
actions are the most relevant spin-spin interactions 
and studying them in a systematic way is necessary in 
order to understand how a local spin excitation would 
propagate through the magnetic media as it is con- 
veyed by each of these interactions. In particular, this 
is relevant in the pump-probe like experiments where 
one is interested in the long-range effect of a (local) 
demagnetization by a laser or any other source of lo- 
cal heating. 

Another example of application of the model con- 
sidered in this work is that of two magnetic layers, 
each represented by its macroscopic magnetic mo- 
ment, separated by a nonmagnetic spacer and coupled 
through the latter by an effective interaction which is 
either EI, DDI, or DMI. The static and dynamic prop- 
erties of such multilayers with effective interactions 
have been studied experimentally and theoretically 
by many groups, see e.g., Refs. 4-6 where ultrathin 
multilayers were studied by the Brillouin Light Scat- 
tering technique and the effective coupling was es- 
timated. This model can also be adapted to a pair 
of macroscopic moments representing two nanopar- 
ticles, embedded in a nonmagnetic matrix. Further- 
more, this somewhat toy problem may serve as a 
benchmark for many-particle systems. Moreover, for 
studying the dynamics of interacting systems one has 
to understand the dynamics of the elementary brick of 
two interacting elements and this can be represented 
as a magnetic dimer (MD). 

One of the questions that we address here is how 
each of these interactions affects the dynamics and, in 
particular, the magnetization reversal of the MD. For 
this we first (semi)-analytically compute the switch- 
ing rates in several situations. Knowing the switch- 
ing time a comparison with e.g. Network Analyzer- 
Ferromagnetic Resonance(NA-FMR) measurements, 
should help us to tell which interaction is most rel- 
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evant in the system studied. 

In the absence of thermal fluctuations, the switch- 
ing is entirely deterministic and occurs at some crit- 
ical value of the effective field, at which the (effec- 
tive) energy barrier vanishes. In the case of DDI the 
switching dynamics has been studied in the past by 
many authors, see for instance Refs.^ ®. For classical 
systems, thermal effects on the dynamics of an MD 
can be accounted for within the Langevin approach 
by solving the system of two coupled Landau-Lifshitz 
equations, one for each magnetic moment, where the 
effective deterministic field is augmented by the (ran- 
dom) Langevin field^. The equivalent approach that 
consists in solving the corresponding Fokker-Planck 
equation (FPE) can also be used^". In the case 
of high-to-intermediate damping regime and high- 
energy barriers, analytical expressions have been ob- 
tained for the exchange-coupled MD^^ using Langer's 
approach^^ and the results have been favorably com- 
pared with the numerical approach based on the solu- 
tion of FPE In Ref it was shown that there exists 
a critical exchange coupling below which the MD re- 
verses its direction via a two-step process, i.e. through 
a fanning mode, whereas above this critical coupling 
the system switches through a coherent mode. At 
the critical coupling, some saddle points become flat 
and Langer's approach, based on the quadratic ex- 
pansion of the energy, ceases to be valid. As a conse- 
quence, the switching rate presents two disconnected 
branches corresponding to the weak and strong cou- 
pling regimes. In Ref the FPE is transformed into 
an eigenvalue problem and the latter is then solved 
using the numerical technique of matrix-continued 
fractions. This work renders a smooth switching rate 
for exchange coupling strengths that compares very 
well with the analytical asymptotes outside the crit- 
ical region. Finally, it is worth mentioning the work 
of Solomon^* where switching processes were studied 
and the longitudinal and transverse switching times 
were computed for a quantum MD at room tempera- 
ture. 

In the present work, we use Langer's approach to 
investigate the dynamics of an MD with three differ- 
ent couplings, EI, DDI, and DMI. The objective here is 
to compare the dynamics and switching modes in each 
case with the aim to provide an answer as to which 
coupling is the most efficient as far as the full switch- 
ing of the MD is concerned and which configuration is 
the most optimal for applications. To accomplish this, 
we consider two orientations of the two anisotropy 
axes with respect to the bond axis, namely the longi- 
tudinal (LA) and transverse (TA) anisotropies. These 
two anisotropy configurations mimic the two cases of 
magnetic films with two limiting thicknesses. In or- 
der to carry the full analytical calculation of the var- 
ious switching rates and provide (approximate) sen- 
sible analytical expressions thereof, the applied mag- 
netic field has been ignored in this work. The other 



reason for this restriction is the need to investigate 
the three couplings and compare the dynamical be- 
havior they entail without the influence of the ap- 
plied magnetic field. In a subsequent work we in- 
tend to include, within a numerical approach, the lat- 
ter and consider more general configurations of the 
anisotropy axes. 

The paper is organized as follows. In the next sec- 
tion we write all the contributions to the MD energy, 
introduce our notation, and present the method for 
computing the various switching rates. In the next 
section, we deal with the three types of interactions. 
We first summarize the previous results for an EI- 
MD and then deal with DDI and DMI, successively. 
We treat the three cases of i) longitudinal anisotropy 
(LA) where both anisotropy axes are parallel to the 
MD bond, ii) transverse anisotropy (TA), with the 
anisotropies perpendicular to the bond direction, and 
iii) mixed anisotropy (MA) with one anisotropy axis 
parallel and the other perpendicular to the bond. In 
each case we investigate the various coupling regimes 
and compute the corresponding switching rates. The 
next section is devoted to a case study of a comparison 
between the three interactions. 



II. ENERGY AND SWITCHING RATE 

In this section we define our notations and provide 
the basic formulae of our calculations, including the 
energy and the switching rate. For the latter we 
briefly summarize Langer's kinetic theory that will 
be used in the case of intermediate-to-high damping 
(IHD) regime. 



A. Energy 

As mentioned in the introduction, the system we 
study consists of two macroscopic magnetic moments 
nii,i = 1,2, which may represent two magnetic lay- 
ers or two magnetic nanoparticles or still two atomic 
moments. Each magnetic moment has an effec- 
tive uniaxial anisotropy that is assumed to result 
from magneto-crystalline and/or shape anisotropy. In 
the case of multi-layered systems, we consider both 
in-plane and out-of-plane anisotropy, thus modeling 
magnetic films with two limiting thicknesses. The 
two layers are coupled via the nonmagnetic layer by 
an effective EI, DDI, or DMI and the corresponding 
dimer will then be referred to as the EI-MD, DDI- 
MD, or DMI-MD, respectively. While exchange cou- 
pling is invariant under global rotation of the system, 
DDI and DMI are anisotropic as they originate from 
a coupling to the lattice. In particular, the magnetic 
state induced by DDI depends on the orientation of 
the vector connecting the magnetic moments, the MD 
bond. For the DDI-MD, for definiteness, we set the 
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latter in the z direction and the corresponding verse 
will be denoted ei2, i.e., ei2 = e^. The connecting 
vector ri2 is then written as ri2 = dei2 where d is 
the distance between the centers of mass of the two 
layers. Since the magnetic layers are assumed to be 
much thinner than the nonmagnetic spacer, the dis- 
tance d is approximately the thickness of the latter. 
The anisotropy axes e^,? = 1, 2 and the applied field 
with the verse e/i are a priory in arbitrary directions 
[see Fig. [l). 




Figure 1: Setup of the DDI-MD with oblique magnetic 
field and arbitrary anisotropy axes. 

In the sequel, we will use spherical coordinates 
for all vectors involved. Hence, for the magnetic 
moments we write = rriiSi, with ||s,;|| = 1 and 
Si{Oi,ipi),i — 1,2. The applied field is written H = 
H Bh, with ||e/j|| 1 and eh{9h, (ph), and the anisotropy 
axes are ei{9^°^\^ ^^t'')- ^ ^ respectively the po- 
lar and azimuthal angles as defined in Fig. [l] Note 
also that the applied field H is to be understood as 
^qH which is counted in Tesla. 

The energy of the MD then reads 



E = Ez+EA + Eint (1) 
where Ez is the Zeeman energy 

= -H • ^ m, = - ^ m, (H • s,) , (2) 

i=l,2 1=1,2 

Ea is the (uniaxial) anisotropy contribution 



(3) 



i=l,2 



The interaction energy _Eint may stem from the ex- 
change (ferromagnetic) coupling 

-Bint = -EexcIi = -^Si • S2, (4) 

from the DM coupling 

Elnt = £'dM = -D • (Si X S2) , (5) 

or from the DDI contribution 

iS,,.. ^ £00, (6) 

with 

2?i2 = 3(^-e,e,-^)-l. (7) 



being the DDI tensor. 

Let's recall that DMI is an anti-symmetrical ex- 
change interaction coming from a combination of low 
symmetry and spin-orbit coupling^^"^^. In the pres- 
ence of disorder, especially at the interface of thin 
films or multilayers, the DMI has been shown to play 
an important role since local symmetry is broken by 
surface effects. Indeed, it leads to large anisotropy 
and may even change the magnetic order, see Ref 
and references therein. In particular, it has been 
shown that DMI is induced by spin-orbit coupling be- 
tween two ferromagnetic layers separated by a para- 
magnetic layer^^. Accordingly, in the present study, it 
is also relevant to investigate its effect on the dynam- 
ics of the MD, on the same footing as the (symmetri- 
cal) effective exchange coupling. 

Investigating the general situation with arbitrary 
orientations for the easy axes is rather involved and 
can only be dealt with numerically. This will be done 
in a subsequent work. In the present work, we choose 
to focus on the qualitative behavior of the various in- 
teractions and investigate how they affect the dynam- 
ics of the system. For this purpose, we consider a sit- 
uation that can be dealt with analytically, thus allow- 
ing for a simpler analysis of the underlying physics. 
More precisely, we assume equal magnitudes for the 
two magnetic moments with equal anisotropies (in di- 
rection and magnitude), i.e. ei || e2 and Ki = K2', no 
external magnetic field. 

In the sequel, we will measure the energy in units 
of the anisotropy energy and thus write 



E 



E 



where 



_ KV 

is the reduced anisotropy energy and also the reduced 
energy barrier in the non-interacting case. Therefore, 
the MD energy reads [the Greek indexes run over 
x, y, z while the Roman indexes run over 1, 2] 



E 

kef 



-2/i^ eh,a ^ Si 



i=l,2 



i=l,2 a, 13 



-Smt (8) 



with 

£lnt — — ^ Sl,Q 



S2,/3- 



(9) 



e"^'' is the fully antisymmetric Levi-Civita tensor of 
rank 3 and is the verse of D. We have introduced 
also the following (dimensionless) parameters 
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_ Hm , _ J 
2KV' " 




(10) 



which imply that all energies are measured in units 
of the anisotropy energy. For instance, h is the usual 
ratio of the magnetic field H to the anisotropy field 

m Ms fia 

B. Relaxation rate 

It was shown in Refs.lSI that Langer's expressionP^ 
for the escape rate from the minimum {9^'"^\ (p^"^'>) 
through the saddle point (^?(^\ takes the more 
compact form 

r=^-^^, (12) 

where Zm and Zg are respectively the partition func- 
tions computed in the vicinity of the minimum and 
the saddle point and \k\ is the attempt frequency 
The latter represents the growth rate of a nucleat- 
ing fluctuation at the saddle point and thus charac- 
terizes the unstable barrier-crossing mode. This ex- 
pression indicates that the escape rate is simply given 
by the ratio of the total current through the saddle 
point to the number of particles (or points in the sys- 
tem phase space) in the metastable state. In fact, 
within Langer's approach the problem of calculating 
the switching rate for a multi-dimensional process is 
reduced to solving a steady-state Fokker-Planck equa- 
tion for the probability density p, i.e., dtp = 0, in 
the immediate neighborhood of the saddle point that 
the system crosses as it goes from a metastable state 
to another state of greater stability. The probability 
density p is connected to the probability current via 
the continuity equation. On the other hand, p can be 
written as peq = e^^"^ /Z times some other function. 
From these two relations one can write the probabil- 
ity current in terms of the partition functioiPSEHI^ gee 
also Ref 21 for great details. Now, since the switch- 
ing rate is given by the total probability flux through 
a surface near the saddle point, Langer's result for 
the escape rate can be achieved by computing the 
energy-Hessian eigenvalues near the saddle points 
and metastable states. From the latter, one then in- 
fers the partition function Zs of the system restricted 
to the region around the saddle point and the parti- 
tion function Z,„ of the region around the metastable 
state. When computing these partition functions, one 



has to identify and take care of each Goldstone mode, 
that is a massless mode or zero-energy fluctuation 
associated with a continuous unbroken global sym- 
metry. The tilde on Zs reminds us of the fact that 
the negative eigenvalue of the energy Hessian corre- 
sponding to the escape route is (formally) taken with 
the absolute value^^. Zs is the product of contribu- 
tions from all eigenvalues. 

In Langer's approach the attempt frequency n is 
computed by linearizing the Landau-Lifshitz equa- 
tion around the saddle point, diagonalizing the result- 
ing transition matrix^^^^, and selecting its negative 
eigenvalue. However, only in a few situations can k 
be obtained analytically. In fact, in the general situ- 
ation, K can only be computed numerically. Accord- 
ingly, one computes the unique^'* negative eigenvalue 
K of the steady-state FPE corresponding to the unsta- 
ble mode at the saddle point as the negative eigen- 
value of the dynamic matrix A/,„„ = — A„(PAfP^)„i„, 
where A„ are the eigenvalues of the energy Hessian at 
the saddle point, M the dynamic matrix, and P is the 
transformation matrix from the initial coordinates to 
the "canonical" ones. 

Therefore, for a given elementary process, i.e., an 
escape from the minimum 

(-^(m)^^(m)) through the 
saddle point (0^"*^ , (f'^''^^), we have to compute the parti- 
tion function 

Z = I (Vs) e-^^(^) 

at the saddle and metastable states. For this, we per- 
form a quadratic expansion of the energy at these sta- 
tionary states. This is where Langer's approach meets 
its limit of validity because such an expansion is only 
meaningful when the stationary point is well defined. 
More precisely, Langer's approach is only valid in the 
case of high energy barriers AE, i.e., when 0AE >• 1 
and also intermediate-to-high dampin^^^IHI. 

In the case of a two-body problem, such as that of 
MD, in the weak coupling regime the magnetization 
of the whole system switches in a two-step process; an 
example is shown in Fig. [2] in the case of EI-MD. The 
first step of switching corresponds to the passage of 
the first magnetic moment from the initial state into 
an intermediate state through the saddle point. This 
step lasts the (switching) time n. The second step 
is taken by the second magnetic moment that then 
proceeds to switch through a second saddle point and 
this step lasts the time t2. The total time required by 
the MD to switch is then t ~ ti + t2, and in terms of 
the switching rate (F = t^^), one has 

11 1 
f " I\ +f^- 

Next, for fully identical magnetic moments one has 
to consider the left-right symmetry and multiply the 
expression above by a factor of 2 leading to the final 
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expression for the switching rate corresponding to the 
two-step process of the MD 



r 



total 



= 2 



(13) 



In fact, one may have other symmetry factors de- 
pending on the system setup. 

Consequently, in the sequel our task will consist in 
analyzing the energy potential surface in each situa- 
tion, studying the various switching paths, and com- 
bining the corresponding switching rates according to 
Eq. ^T3\ . For each elementary step corresponding to 
an escape from a minimum through a saddle point we 
will use Eq. ( [T2| to compute the corresponding switch- 
ing rate. 

Defining the characteristic time of the underlying 
material ts = {-/HAy^ = Ha/C^^Ka), where 7 ~ 
1.76 X 10^^ (T.s)^^ is the gyromagnetic factor, the final 
(dimensionless) switching rate may be given in s^^ 
upon multiplying by 



ikeT ^ 1 / 2-/KV \ f kBT\ 
Ma 2 Ms y \KV) 



2 cr ' 



(14) 



For cobalt, for instance, we have fia — 1-57 x 
10-^^ Am-'^atom-'^, Ka = 2.53 x 10-^'^Jatom-\ lead- 
ing to ts = 1.76 X 10-".s. 

In this work we compute the relaxation rate by com- 
bining Langer's approach, which is valid in the IHD 
regime, and the Landau-Lifshitz equation (LLE) with 
its phenomenological damping parameter. Using the 
LLE for obtaining the attempt frequency in the pref- 
actor in Eq. ( [12| leads to spurious effects when one 
formally takes the limits a — > or a — >■ ooW How- 
ever, Langer's expression for the relaxation rate may 
be "regularised" by using Gilbert's damping instead 
of the Landau-Lifshitz damping. Indeed, the latter 
may be shown to be identical to Gilbert's if the gy- 
romagnetic factor 7 is replaced by 7* = 7/(1 + 0^). 
In the present calculations this amounts to replacing 
the scaling time tg defined above by t* = (l + a^). 
In fact, Landau-Lifshitz and Gilbert's equations are 
related by the transformation 



III. EXCHANGE COUPLED MAGNETIC DIMER 

In this section, for later use, we briefly summarize 
the results of Ref.'^. The situation is sketched in Fig. 

m 

In the case of parallel easy axes and longitudinal 
field, it was found that there is a critical exchange 
coupling jc that depends on the applied field and 
anisotropy constant, i.e. jc{H,K), above which the 
MD behaves as a macrospin with a double energy bar- 
rier that switches from a metastable state to a more 
stable one in a coherent manner, see Fig. |2](right). Be- 
low jc the system is weakly coupled and switches in a 
two-step process through two different escape routes 
(saddle points), see Fig [2 (left). Setting eh \\ ei || 62, 
in the notations of Eq. ( 10 1, the critical exchange cou- 
pling jc was found to be 



Jc = 1 - 



(15) 



For j > jc, the MD switches from the metastable 
state (tt, tt, (p), that is a ferromagnet, against the field, 
into the ferromagnetic state (0, 0, </?), through the sad- 
dle point (arccos {~h) , arccos {—h) , tp). The angle tp is 
arbitrary because of the uniaxial symmetry and the 
equality of the spin polar angles is due to the fact 
that the two spins are identical (same amplitude and 
same anisotropy). The corresponding switching rate 
is given by (in zero field) 



2^ 1 + 

7T v/r^v7 



(16) 



In the weak-coupling regime j < jc, the switching 
rate is obtained by combining the switching rates cor- 
responding to the two escape routes taken by the two 
spins, see Fig. |2](left). The result is somewhat more 
involved and given in Refs.^^, see also RefP. 

The escape rate for the El-coupled MD will be com- 
pared to the other two cases of DDI- and DMI-coupled 
MD. For this purpose, we recall here the energy bar- 
riers (in the absence of the magnetic field) 



(17) 



7^ 



1 



1 + 



We recall that our main objectives in this work are: 
i) an investigation of the behaviour of the relaxation 
rate as a function of the MD coupling and ii) a pair- 
wise comparison of three types of layer coupling (ex- 
change, dipolar, and Dzyaloshiski-Moriya). In partic- 
ular, we do not investigate the damping dependence of 
the relaxation rate. A thorough study of all crossovers 
between the various damping regimes is given in Ref 
|20 where one can see that the boundaries between the 
regimes are not simply a — 1. 



IV, DIPOLAR COUPLED MAGNETIC DIMER 

We consider both the longitudinal and transverse 
anisotropies, i.e. ei || 62 || ei2 || and ei || 62 -L ei2 | 
e^, respectively, which will be referred to as the LA 
and TA setup, respectively. Furthermore, we will also 
discuss the case of mixed anisotropy, ei || 812, e2 -L 812, 
which would mimic the case of an MD with a suffi- 
ciently thin film coupled to a sufficiently thick one. 
Comparison with the exchange-coupled MD will be 
done only in the case of LA, considered in Refs.^^. 
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A few general assumptions allow us to simplify the 
problem without any loss of generality as far as the 
underlying physics is concerned. Indeed, in the se- 
quel, we will assume the following. If there is no mag- 
netic field, the equilibrium orientation of the net mag- 
netic moment is in the plane defined by the DDI axis 
and the anisotropy axes. The latter are assumed to 
lie in the xz plane, i.e., ip'f^ ~ ^^^^ = 0- Consequently, 
the energy becomes 



-a ^ cos^ 



g(a) 



1=1,2 



(18) 



— cr^ [2 cos 6*1 cos d-i — sin sin 6*2 cos ((ySi — Lp2)\. 

The analytical study will be further restricted to 
the following three cases: 

1. Longitudinal anisotropy (LA): both anisotropy 
axes are parallel to the MD axis ei2, i.e. (f-f' — 

= Q^2^ . Moreover, due to the fact that all con- 
tributing fields are acting in the same plane, the 
two magnetic moments of the MD will move in 
the same plane so that (^i — 

2. Transverse anisotropy (TA): both anisotropy 
axes are perpendicular to the MD axis, 6*^"' — 

TT _ n(a) 

1 —^2 ■ 

3. Mixed anisotropy (MA): one anisotropy axis is 
parallel to the MD axis and the other perpendic- 
ular to it, 61^"' = 0,6*^"^ = f. 

Differentiating with respect to the remaining vari- 
ables, i.e. the two polar angles O^.i — 1,2, leads to 



the various equations for the stationary states whose 
solutions depend on the anisotropy setup. 



A. Longitudinal anisotropy 

The whole set of stationary states is given by 

(01,^2) = (0,±7r),(±7r,0),(±^,±7r),(±7r,T7r), 

(±2'±2)'K'^2) 

together with the following ones (which are saddle 
points) for ^ < I or ^ > 2, 

(cos0i,cos02) = (±a;^,±a;2) , (^xX^^xr^) (19) 

where 

1 L 3 




It can be checked that the radicant is always pos- 
itive and no additional special ranges are found re- 
garding the existence of the roots. However, as x\ 
and x\ are cosines they must satisfy — 1 < a;- < 1, 
which is only true for ^ < |. Hence we can identify 
two different regimes, the weak-coupling (WC) regime 
^ < |, and the strong-coupling (SC) regime C > §■ We 
will see later that this critical value corresponds to 
the vanishing of the smallest eigenvalue of the energy 
Hessian at one of the energy minima. It also marks 
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(c) 5 = 0.6 (d) ^ = 1 

Figure 3: Evolution of the energy potential surface of a DDI-MD with longitudinal anisotropy configuration as 

^ increases, with a = 1.5. 



the nucleation of a particular switching mode and al- 
lows us to determine the nucleation field^ ^. 

The energy potential surface for this situation is 
shown in Fig. |3]for a varying DDI strength ^. 

1. Weak coupling (£, < p 

The minima correspond to ferromagnetic (FM) 
states along the DDI axis 

(01,02) = (0,0), (tt,^). (21) 

The metastable states are the anti-ferromagnetic 
(AFM) states along the DDI axis 

(01,02) = (0,^),(^,0) (22) 

and the maxima are the (anti)ferromagnetic states 
perpendicular to the DDI axis 

= (±i,±|),(±i.^i). (23. 



Finally, the saddle points are located at 

(^1,^2) — (ei arccos (£22^1) 7 £1 arccos (£22:2)) (24) 

where the signs £1 = ± and £2 = T are independent 
of each other. 

Now we compute the switching rate in this cou- 
pling regime. In Fig. [3^ we see that the system goes 
through the following steps: i) from the state (tt, tt) to 
the state (0, tt) through the saddle point 

(cos0i^^ = a;+,cos0^^^ = -a;+) (25) 

[see Eq. ([24]l], and then ii) it passes from the state 
(0, vr) to the state (0, 0) through the saddle point 

(^cos0f^ = a;j;,cos0*^^ = -x-^ . (26) 

These transitions are sketched in Fig. |4](left). 

Therefore, in order to compute the switching 
rate corresponding to the two-step process (tt, tt) 
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(0, vr) (0, 0) we need to compute the switching rate 
of each step and combine them according to the rule 
in Eq. (|13[l where the individual switching rates are 



then computed using Langer's expression ( 12 1. 

In order to compute the switching rate Tr)^(o.7r) 
we first compute the partition function at the mini- 
mum (tt, tt) and at the saddle point (9[^^ , O^^'^ ) , namely 



Likewise, the partition function at the metastable 
state (0, tt) is 



(2^) 



=2<t(1-C) 



a2 (2 - e) (2 - 30 



(28) 



This result can also be found upon noting that since 
the anisotropy is uniaxial, in order to change the 
system minimum from the ferromagnetic state (tt, tt) 
to the anti-ferromagnetic state (0,7r), we can simply 
change the sign of the interaction, i.e. replace ^ by 

The lowest eigenvalue of the energy Hessian at the 
metastable minimum is A = cr (2 — 3^. As mentioned 
earlier, we see that when this eigenvalue vanishes it 
yields the critical value for the DDI coupling, namely 
^ = 2/3. Indeed, the nucleation field in this case is 
hn oc CT (2 — 3^ which coincides with the result of Refs. 
I7I8I10I In Ref \W\ the authors defined the parameter 
for the DDI strength fcjnt as the ratio of the DDI field 
to the anisotropy field. In our notations, feint = C/2- 
The critical value of feint — 1/3 coincides with our con- 
dition C 2/3. 

At the saddle point (|25| the energy is 



(1) 



(29) 



and the energy barrier separating the minimum (tt, tt) 
from the saddle point ^?2^'^ is 



(30) 



This is plotted in Fig. [5] where we see that the 
(weak) DDI brings a correction a (2^ + |^^) to the free 
MD energy barrier Ae ~ a. For ferromagnetic order 
(^ > 0), this correction enhances the energy barrier 
and this is compatible with the fact that due to the 
ferromagnetic coupling, it is more difficult for the first 
spin to switch. 

Within Langer's approach to the calculation of the 
switching rate the saddle point may retain a subgroup 
of the symmetry group of the system, in which case 
some of the Hessian eigenvalues at the saddle point 
vanish, and then a special treatment is required for 



this situation. For the one-spin problem with uni- 
axial anisotropy, for instance, the saddle point has 
a C/(l) symmetry around the z axis, or with respect 
to the rotation TZ{e2,ip). This leads to a vanishing 
eigenvalue of the Hessian at the saddle point corre- 
sponding to fluctuations with respect to the angle ip 
(~ if^ + p), where ip^ is the value of the azimuthal an- 
gle at the saddle point. Likewise, the saddle points 
( p5| [26| l have rotational symmetry with respect to the 
azimuthal angle and as such one should use the en- 
ergy 

£ = -a [cos^ 6i + cos^ 62] 

—a£, [2 cos 9i cos 62 — sin di sin 82 cos ip] 

where ip = (pi — (^2- The saddle points ( (25|[26| should 
refer both to the polar angles 6i and to a given value 
of the azimuthal angle tp, even if it is arbitrary. 
Let us then perform the expansion around the sad- 



die point (j) 



and expand the energy above to second order in the 
small variables ti,pi. Doing so we obtain 



a 

-a 



where 7r± = 



+ 2 

Pi ± P2 





e 



Note that the mode corresponding to 7r+ is a soft 
mode, i.e., a zero-energy mode, or still a Goldstone 
mode. This corresponds to the U (l) symmetry men- 
tioned earlier and that must be dealt with properly in 
order to avoid the divergence of the partition function. 
In fact, integration over this variable simply yields 
the factor 2n. The partition function is then given by 



-e "0 



(31) 



Hence, gathering the results in Eqs. < \27\ and (311 
into the first equation of p2]l, we obtain 



(7r,Tr)^(0,Tr) 



.(1) 



a (2 + 30(2 + 0^_A£(i' 



2a 



At the second saddle point ( [26| l the energy is the 
same, i.e.. 



c(2) _ 
^0 — 
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Figure 4: Escape route of the MD with LA. 



whereas the energy barrier separating the minimum 

(2) (2) 

(0, vr) from the saddle point {6\ ',62 ) reads 



a 1 - 2e + -r 



(32) 



Here the correction to the energy barrier of the free 
MD starts with a negative coefficient which imphes 
that this energy barrier decreases when ^ increases, 
see Fig |5] Indeed, as the DDI becomes stronger, the 
first moment, which has already switched, exerts a 
stronger force on the first moment. On the opposite. 



the second energy barrier (30 1 increases with ^. When 



the system is in the initial ferromagnetic state (tt, vr), 
it is much more difficult for the first magnetic moment 
to break free from the ferromagnetic coupling. 
The corresponding switching rate then reads 



(0,Tr)^(0,0) 



.(2) 



a (2 - 30 (2 



2a 



Note that we have the symmetry F 



(0,7r)^(0,0) 



r(7r,7r)-i.(o,7r) — > — 0- I^i both switchiug rates the 
prefactor k has been computed numerically. The to- 
tal switching rate of the MD is obtained upon using 
Eq. ([13| and is plotted in Fig. (6] 



2. Strong coupling ((.>2/3) 

The situation now is sketched in Fig. |4] (right) with 
the minima 



71, P2 



) = (0,0),(7r,^) 



maxima 



(33) 



(34) 



^)=(±2'^2 



(35) 



Comparing with the WC regime, we note that the 
global minima remain the same, while the metastable 
minima and the saddle points are no longer station- 
ary states. On the other hand, the anti-ferromagnetic 
states given by (±|, =f|) change as well, being max- 
ima for the WC regime they turn into saddle points 
in the SC regime, creating a switching path in which 
both magnetic moments switch coherently in a one- 
step process, as shown in Fig. |4] 

Following the same procedure as for WC, we ob- 
tain the expression for the switching rate in this SC 
regime ofDDI-MD 



Tlasc 
where 

K = ( 1 - 



aK\l — (3^ 




2 + ^ 



--(2+0 (36) 




and saddle points 



(37) 



Fig. [5] shows the evolution of the energy barrier as 
a function of the DDI coupling ^. We clearly see the 
two coupling regimes separated by the critical value 
^ = 2/3. In the WC regime we do see the two different 
energy barriers corresponding to the two steps of the 
reversal process. The energy barrier is continuous for 
all ^, including the critical region. Moreover, the sec- 
ond step disappears at ^ = 2/3 and the dynamics of 
the system becomes a one-step process. 

Figure |6] shows the behavior of the reduced switch- 
ing time T — 1/T — t/2ts as a function of the reduced 
anisotropy barrier a for both coupling regimes. 

The total switching time may be written as 
Twc = ■^(7r,7r)^(o,7r) + T(o,7r)^(o,o) sud it is clear that 
T(n,n)^{o,n) > T(o,7r)^(o,o) siuce iu the first step one 
of the magnetic moments has to win against the ef- 
fective FM coupling. Consequently, as the DDI cou- 
pling increases the magnetic moment switching dur- 
ing the first step becomes more and more difficult to 
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Figure 5: Energy barrier as function of ^ for the 
LA-MD, with = 3/2. 




2 4 6 8 10 

KV/kgl 

Figure 6: Reduced switching time r of the DDI-MD 
with LA as a function of a for weak and strong 
couphng regimes . 

achieve and thereby the corresponding switching time 
increases, which explains why tlawc < tlasc- 

B. Transverse anisotropy 
In spherical coordinates the energy now reads 



£ = -a (sin^ 6ii + sin^ (38) 
— (7^ [cos (6'i + Q-i) + cos 6*1 cos ^2] 

In this situation, we find that there are two minima 
with the corresponding lowest eigenvalues cr(2 — ^) 
and cr (2 — 3^) whose vanishing leads to the two criti- 
cal DDI couplings ^ = 2 and C = 2/3. 

1. Weak coupling^ < 2/3) 

The absolute minima of the system are now 

(01,02) = (±|,t|) (39) 
while the local minima are given by 

(01,02) =(±|,±|). (40) 

On the other hand, the absolute maxima are the 
anti-ferromagnetic states 

(01,02) = (0,^),(^,0) (41) 

while the local maxima correspond to the ferromag- 
netic states 

(01,02) = (0,0), (tt,^). (42) 
Finally, the saddle points are 

(^'1,^2) = (ei arccos(e22;i), ei arccos(£2a;2)) • (43) 

The escape route is sketched in Fig. [7] (left). 

In the present case, as can be seen in Fig. [8] 
the initial state is given by (|,— |), the metastable 
minimum is (f , f ), and the first and second sad- 
dle points are given by (arccossf ,arccosa:;J) and 
(arccosX]^, aiccosa;2^), respectively. In this weak cou- 
pling regime, the MD switching is again a two-step 
process and the corresponding switching rates are 
given by 



-(1) 

TAWC 

.(2) 
TAWC 



TAWC 



2^/37l 



V(3e + 2)(l + e) (2-0 X e-¥(f+«)(i-|) 



(44) 



where the attempt frequencies k^*', i = 1,2 are calcu- 



lated numerically; a is defined in Eq. ( |20| l. The total 
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Figure 7 : Escape route for the three DDI coupling regimes with TA configuration 




switching rate is given by Eq. p^ . 

2. Medium coupling f| < ^ < 2J 
The minima are given by 



, TT vr 

71,(72) = ±-,T77 



(45) 



the maxima correspond to the anti-ferromagnetic 
states 



(01,02) = (0,^),(^,0), (46) 
and the saddle points are the ferromagnetic states 



and 



= (0,0), (tt,^) 



(47) 



(48) 



Starting from the initial state (|,— |), the switch- 
ing is a single-step process in which the two magnetic 
moments switch coherently through the saddle point 
(0, 0) leading into the state (— |, |), as shown in Fig. 
|7] (middle). 

The switching rate for this coupling regime reads 



TAMC 



«(C-2) / (3^ + 2)(l + -.(2-e) 
TT V (3C- 2) 



(49) 



where the attempt frequency has been obtained ana- 
lytically and is given by k = 2q: (^ - 2). 



3. Strong coupling (£, > 2) 

Here the DDI field is twice larger than the 
anisotropy field and the minima are the ferromag- 
netic states 

(01,02) = (0,0), (tt,^), (50) 
the maxima are the anti-ferromagnetic states 

(01,02) = (0,^),(7r,0), (51) 
and the saddle points are located at 

(..».)^(4±f).(±f.^f)^ 152, 

As the DDI coupling increases and the system en- 
ters the strong coupling regime, the states that pre- 
viously were minima in the medium coupling regime 
become saddle points, and vice versa. The system still 
relaxes in a one step process, but this time from the 
initial state (0,0) through the saddle point (— f,|). 
This is sketched in Fig JT] (right). The switching rate 
in this case reads 



TASC 



(3^ - 2) 

(2 + 30(1 + 0' 



= -'^(«-2) 



(53) 



with 
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Figure 9: Energy barrier of the MD with TA as a 
function of ^, cr = 3. 



Figure 10: Reduced switching time r of the MD with 
TA for the three different couphng regimes as a 
function of a. 



more complex behavior, as can be seen in Fig. 11 



Fig. [9]shows the evolution of the energy barrier as ^ 
changes. Similarly to Fig. |5] the different regimes are 
clearly identified, so are the continuity of the barrier 
and the disappearance of the second step at ^ = 2/3. 
At ^ = 2, we can see that the energy barrier vanishes 
forming a "furrow" that connects directly (with no en- 
ergy barriers) the states (0,0), (tTjTt), and (±|,=Ff)- 
Figure irol shows the behavior of the switching time in 



the three regimes of DDI coupling, namely WC, MC 
and SC; it is a plot of the inverse of the expressions 
( |44| , ( |49| , and ( (53) l, respectively. 

We again observe that the SC switching time is 
longer than that of the WC, similar to what we have 
observed in the case of LA, see Figj6] Furthermore, we 
see that, as a function of cr ^ 1 /T, the (logarithm of) 
the switching time in the WC regime is not a straight 
line. This implies that the prefactor plays a dominant 
role. As ^ increases the switching time becomes dom- 
inated by the Arrhenius (exponential) law where the 
prefactor is a constant thus leading to a straight line 
in a logarithmic plot of the switching rate as a func- 
tion of cr. One of the consequences of a dominant pref- 
actor is that the switching rate becomes quite sensi- 
tive to damping and thus to the coupling of the system 
to its thermal bath and to the various fluctuations. 



Analysis of the stationary states reveals the presence 
of two coupling regimes as well, but with a critical 
value now at ^ = As before, switching in the WC 
regime is again a two-step process while it is a one- 
step process in the SC regime. The stationary states 
are given by Eq. |Al| together with 

{01,92) = (0,±7r),(±7r,0),(±^,±7r), 

(±-,T.),(±2'±2)'(±2'^2)- 

It is worth investigating this case because it corre- 
sponds to another interesting situation with a rather 
thick layer coupled to a rather thin layer. It may also 
be relevant in the situation of a soft magnetic layer 
coupled to another hard magnetic layer. 

The remaining stationary states are rather cumber- 
some and they are relegated to the appendix. 



1. Weak coupling((, < 

The individual switching rates are given by Eq. 
(A2]l, and the full switching rate is given by Eq. 
( 13|l. The corresponding analytical expressions are 
too cumbersome and are thus given in the appendix. 



They are plotted in Fig. 13 In this case too, k is com- 
puted numerically. 



C. Mixed anisotropy 

2. Strong coupling > 

Whereas the LA and TA configurations for DDI 
coupling have a somewhat similar behavior in terms The minima now are the FM states 
of energy potential surface evolution and coupling 

regimes, the MA shows a completely different and (6*1,^2) = (0, 0) , (tt, tt) 



13 




(a) Weak Coupling 5 = 0.2 (b) Approximate escape route (c) Strong Coupling g = 1-66 

WC 

Figure 11: Energy potential surface for the different regimes of the MD and WC escape route with MA and 



and the maxima are the AFM states 

(01,02) -(0,^), (71,0) 

while the saddle points are 

Starting in the WC regime, as the DDI coupling in- 
creases, the minima given by 

(±(arccosa;^ — arccosxI)/2,±(arccosa;;jL + arccos a:;I)/2) 

start to merge into the saddle point (0,0). When the 
system enters the SC coupling regime, these two min- 
ima completely merge leading to the transformation 
of the state (0, 0) from a saddle point into a minimum. 
The minima 




Figure 12: Changes of the energy barrier as a 
function of ^ with MA and a — i. 



(±7r =F (arccos + arccos a;_)/2, 
±7r =F (arccos — arccos a;I)/2) 

have a similar behavior around the state (tt, vr), while 
the saddle points (±f,=Ff) do not change. The es- 
cape route in the SC regime is similar to that of the 
SC for LA and TA. However, the shape of the energy 
potential surface presents curved paths instead of the 
usual straight paths, see Fig. [Tl] 

The escape rate from the initial state (0, 0) through 
the saddle point (— f , |) is given by 



MASC 





(l + 2e)'-r2 


[4^2 _ rl] 






[(C + 2r)|^-2r|] 



,(54) 



where r = vTTC^, Vp = VTT^. n is computed nu- 
merically. 



t/2ts{.a = 5) 


LA 


TA 


MA 


WC 


355.83 


1.48 


20.29 


SC 


1.04 X 10*^ 


11.57 


505.99 



Table I: Reduced switching times for the WC and SC 

for the three anisotropy configurations. ^ = 0.3 for 
WC. e = 1.5 for LA-SC and MA-SC. ^ = 2.5 for TA-SC 



Fig. [12] shows the energy barrier as a function of ^. 
The energy barrier along the second step of the pro- 
cess disappears at the critical value C = 2/ \/3, where 
the dynamics of the system becomes a one-step pro- 
cess. As the DDI increases to higher values, the en- 
ergy barrier increases, leading to a constant increase 
in the switching time for high ^, as can be seen in Fig. 



13 where the switching time for SC is much higher 
than that ofWC. 

By way of comparison we gather in Table |I] some 
results for the DDI-MD. From this table, for cr = 5, we 
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Figure 13: Reduced switching time t of the DDI-MD 
with MA for weak and strong coupling regimes. 



t(a = 8, C = 0.346) 


LA 


TA 


MA 


3.42 X IQ-^s 


9.84 X W~^°s 


7.07 X IQ-^s 



Table II: Switching times for the atomic Cobalt 
dimer for each of the three anisotropy configurations. 



see that for the three anisotropy configurations the 
switching time increases then decreases with the DDI 
coupling ^. The fastest switching seems to occur for 
weak coupling and transverse anisotropy. This means 
that a magnetic dimer composed of two rather thin 
films should exhibit the fastest dynamics. 

Table [n] shows the switching times for cobalt atoms 
in the three different anisotropies. The param- 
eters used for the calculations are fia — 1-57 x 
10-^^ Am-^atom-^, Ka = 2.53 x 10-'^^Jatom-\ d = 
2 X 1.52 X 10-i°m and = 1.76 x 10""s. 

It is clearly seen that in this case TA leads to the 
shortest switching time. As such, in a chain of atoms 
a magnetic excitation should propagate faster if the 
anisotropy is normal to the chain axis. 



V. DM COUPLED MAGNETIC DIMER 



As discussed in section II A DMI is also relevant 



in the present study and is investigated on the same 
footing as EI and DDI. Its effect is compared to that 
of the latter on the switching mechanisms of the MD. 
In order to investigate the effect of pure DMI, we con- 
sider the energy in Eq. ^ without the magnetic field 
and without the EI and DDI. In Ref it was shown 
that for a simple cubic lattice, on the (10 0) surface 
the DMI vector D lies in the layer plane and thus in- 
duces perpendicular anisotropy. Accordingly, in Eq. 
^ we drop the Zeeman energy, the EI and DDI con- 
tributions. We consider two situations where the DMI 



vector D lies in the MD plane and the anisotropy easy 
axes parallel or perpendicular to it. 



A. D parallel to the anisotropy axes 

After simplification, the reduced energy is a func- 
tion of only the polar angles 9i , 62 and is given by 

£{^1,02) = -fc(cos^0i +cos^6'2) + (5 sin 6*1 sin 6*2. (55) 

The anisotropy parameter k = 0, 1 is simply a 
switch introduced so as to be able to keep track of 
the anisotropy contribution in the subsequent results. 
Analysis of the stationary points yields 



e(0,0) = £(7r,7r) = e(0, tt) = e(7r, 0) 



e(± — , =F arcsin 



) = -* 



From this analysis, we find that there is a critical 
value for the DMI which separates the weak and the 
strong coupling regimes. In our normalization with 
respect to the anisotropy energy [see Eq. ^T0\ ] this 
critical value is 6/k ~ 2, see Fig. 15 (left). 



1. Weak coupling S/k < 2 

The minima are the FM and AFM states in the di- 
rection of anisotropy 

(01,02)- (0,0), (^,7r),(0,7r),(^,0), 

the maxima are at 

and the saddle points are 



2k 

2fc / '^2 



The energy potential surface is shown in Fig. [14] 
(a). It can be seen that the net magnetic moment goes 
from, say the minimum (0,0) to the minimum (tTjTt) 
through a two-step process that can proceed along two 
symmetrical paths. Each one of these goes through 
the first saddle point (—§, arcsin (^)), passes into 
the local minimum (tt, O) and crosses the saddle point 
(-^ -arcsin (4) ,-f). 

The switching rates corresponding to these two 
steps have the expression 



^f(4-5^) 



i = 1,2. 



(56) 
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(a) Weak Coupling 5 = 1 (b) Strong Coupling 5 = 3 

Figure 14: Energy potential surface of the different regimes for the DMI-MD with LA. 



where the attempt frequencies are computed nu- 
merically. Upon counting the symmetry factors, we 
obtain the total switching rate 



j^DMI _ o I 
LAWC — ^ I 



(57) 



2. Strong Coupling 

In this regime, the DMI wins against the anisotropy 
field leading to a minimum with perpendicular mag- 
netic moments, lying in the plane normal to the 
anisotropy axes since the DMI vector D is oriented 
along the latter. As such, the minima are 



, vr TT 
71,^2) = ±77'T 



"2' 2 



while the states 



are maxima and the saddle points now are 

(01,02) = (0,0),(7r,^),(0,^) (7r,0). 



(58) 



Hence, the system may escape from the state 
(— |, f) into the state (f f), thus reversing its re- 
sultant magnetic moment, along two different paths 
comprising the saddle points (0,0) and (7r,7r). The 
switching rate of escape via one of these paths is given 
by 



T^DMI 
^ LASC 




-cr(5-2) 



(59) 



As is usually the case in the SC regime, the attempt 
frequency has been obtained analytically. 

Fig. 14 shows the 3-D energy potential surfaces for 
both WC and SC regimes. In Fig. [15] we plot the 



energy barrier (left) and the switching time (right) 
for the DMI-MD in the case of D parallel to the 
anisotropy axes. It is seen that the energy barriers 
for the two steps in the WC regime are equal and 
decrease quadratically with the DMI strength S [see 
the argument of the exponential in Eq. (|57|l]. At the 
critical value of the DMI coupling, 5c/ k = 2, the en- 
ergy barrier vanishes and immediately after that it 
increases linearly with 5/k, as can be seen in Eq. ( 59 1. 

Here again we see that the stronger is the coupling 
the slower is the MD switching. Note that in this 
regime, the saddle point corresponds to the state with 
the two magnetic moments along the easy axes and, 
more importantly, parallel to the DM vector D. To go 
through this saddle point the magnetic moments have 
to break free from the interaction and also to circum- 
vent the anisotropy due to the DM interaction. This 
implies that the DMI leads to a longer switching time 
than the EI [see below]. Indeed, in the latter case 
switching is achieved against the (exchange) coupling 
while in the former it is achieved against the (DMI) 
coupling and the induced anisotropy. 



VI. MOST EFFICIENT COUPLING IN A 
MAGNETIC DIMER 

In this section we present a pairwise comparison 
of the different interactions with regard to their ef- 
fect on the MD switching and on the corresponding 
switching time. On one hand, we have the short- 
range interactions EI and DMI, symmetric and anti- 
symmetric, respectively. On the other hand, we have 
the long-range and antisymmetric interaction DDI. 
We first compare the EI with DDI and investigate 
the effects pertaining to the MD bond. Then, comes 
the comparison between the spins scalar-product and 
vector-product interactions, i.e. EI and DMI. Finally, 
we compare the DDI and DMI. 
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Figure 15: Energy Barrier (left) and reduced switching time (right) of the DMI magnetic dimer with LA for 

WC and SC. 



A. EI versus DDI 

For the EI-MD, when the exchange coupHng ex- 
ceeds the critical value the energy barrier [see Fig. 
2 of Ref ^^J becomes independent of the exchange cou- 
pling as soon as the saturated ferromagnetic state is 
reached. For the DDI-MD, the situation is fundamen- 
tally different because the energy barrier continues 
to increase as the DDI coupling increases, see e.g. Eq. 
( [36| for the SC regime, where A£ = + This 
isdue to the fact that the distance between the two 
magnetic moments, belonging to the layers or to the 
magnetic nanop articles, plays a crucial role. Indeed, 
this distance cannot be smaller than a certain min- 
imal value that corresponds to the thickness of the 
nonmagnetic spacer (in the case of two magnetic lay- 
ers), or to the sum of radii of the two particles, or to 
the inter-atomic distance. Therefore, it is understood 
that the ^ axis must be cut off at a given value be- 
cause the unlimited increase of f simply reflects the 
unphysical asymptotic limit d — > 0. 

Let us now compare the switching times of the MD 
with LA when coupled via EI or DDI, in both the WC 



and SC regimes. The results are shown in Fig. 16 
where the (reduced) switching time is plotted as a 
function of cr = KV/k^T. For these calculations, both 
the EI-MD and DDI-MD switch from the same initial 
state (0, 0) into the same final state (tt, tt). 

As was discussed earlier, apart from the fact that 
the switching time obviously increases with a (or with 
decreasing temperature) for both EI-MD and DDI- 
MD, we see that for both coupling regimes there is 
a critical value (Tc at which the switching times corre- 
sponding to EI and DDI intersect each other. Indeed, 
the DDI is always faster than the EI for low values of 
(T (below (Tc) and the situation reverses for values of 
a higher than (7^. The EI energy barrier is constant 
while that of DDI continues to grow. So, below dc the 



prefactor of the switching time prevails and the DDI 
is more favorable for a fast switching. However, as 
CTc is exceeded, the energy barrier prevails over the 
prefactor and thereby the ever growing DDI energy 
barrier leads to a slower switching than via EI. The 
expression of ctc is obtained in terms of the ratio of 
the switching rates in Eqs. and ( (36) l, which is of 
the form Prefactor/e'^^ with 



Prefactor 




where 



-l)(e + 2)3g + 2 
^ (3^ - 2) j + l 



More precisely, we have 



e + 1 



In (Prefactor) 



The critical value CTc is a decreasing function of the 
ratio ^/j, which is simply due to the fact that the 
stronger the DDI the smaller is at which the en- 
ergy barrier prevails over the prefactor of the switch- 
ing time. 

In conclusion, at low temperature, the EI-MD 
switches faster than DDI-MD. 



B. EI versus DMI 



Fig. 17 shows a comparison between the reduced 



switching times of the EI-MD and the DMI-MD in the 
weak coupling regime. 
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Figure 16: Reduced switching time, r, versus a ~ KV/kBT for the EI- and DDI-MD, in the absence of the 
magnetic field, for weak-couphng regime (left) and strong-coupling regime (right). 





Figure 17: Reduced switching time, r, versus 
CT = KV/kBT for the EI and DMI-MD, in the absence 
of the magnetic field, for weak-coupling regime. 



Figure 18: Reduced switching time, r, versus 
a = KV/ksT for the DDI-MD and DMI-MD, in the 
absence of the magnetic field, for weak-coupling 
regime. 



The initial and final states are identical for both in- 
teractions and are given by {01,62) ~ (tTjTt) (initial 
state), (6*1, 6*2) = (0,0) (final state). 

Unlike the previous cases (Fig. |16| l there is no in- 
tersection of the two curves at finite values of ct in a 
range where our approach is applicable, i.e. where the 
second-order expansion of the energy is applicable. 

From these results we see that the EI-MD has a 
shorter switching time than the DMI-MD. Indeed, in 
the weak coupling regime, for the dimer to switch one 
of its magnetic moments has to cross a saddle point 
into an intermediate state. In so doing, it has to cir- 
cumvent the energy barrier associated with its cou- 
pling to the second moment. In the case of the DM 
coupling, in addition to the breaking free from the 
coupling there is a constraint related with the orien- 
tation imposed by the DM vector D, i.e. the inherent 



anisotropy. 



C. DDI versus DMI 



In Fig. 18 we compare the reduced switching times 



of the DDI-MD and the DMI-MD in the weak coupling 
regime. Two different curves for each interaction are 
presented. 

The initial and final state are identical for both in- 
teractions and are given by (^1,^2) = {t^,t^) (initial 
state), [61,92) = (0,0) (final state). 

We first see that there exists a Oc and that DDI 
leads to a faster switching than the DMI for a < a^- 

Switching in both cases is performed against the 
spin coupling and the anisotropy. However, the latter 
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has a stronger effect in the case of the DMI dimer. 



VII. CONCLUSION 

We have considered a magnetic dimer as a model 
system of two magnetic moments, atomic or macro- 
scopic, coupled by either exchange, dipole-dipole, or 
Dzyalozhinski-Moriya interaction. This is a quite 
general system since the two magnetic moments may 
be those of two thin layers coupled by an effective in- 
teraction through a non magnetic spacer, two mag- 
netic nanoparticles in a hosting matrix or on a sub- 
strate, or still two atomic magnetic moments. We 
have identified various coupling regimes and investi- 
gated the switching mechanisms of the system in each 
regime and in different anisotropy configurations. In 
each situation, we have computed the energy barrier 
and, for the high-to-intermediate damping, we used 
Langer's approach to compute the switching rate and, 
in some cases, provided the corresponding analytical 
expressions. 

We have investigated how the energy barriers are 
affected by the coupling. For instance, for the dipole- 
dipole interaction we find that the energy barrier may 
either increase or decrease with the coupling depend- 
ing on the coupling regime. In the weak-coupling 
regime, we find that the switching rate, as a func- 
tion of temperature, does not follow the simple Ar- 
rhenius law because the prefactor dominates over the 
exponential. Furthermore, transverse anisotropy or 
equivalently, rather thin magnetic films, seem to ex- 
hibit the fastest switching process, as compared with 
the longitudinal or mixed anisotropy. 

We then compared the three interactions with 
regard to their efficiency in switching the mag- 
netic dimer. Comparing exchange and DDI led to 
the conclusion that below some critical temperature 
the exchange-coupled MD switches faster than the 
dipolar-coupled MD. Next, comparing the isotropic 
and anisotropic exchange and Dzyaloshinskii-Moriya 
interactions we have seen that in the latter case the 
inherent anisotropy makes the switching longer. Al- 
together, we have 

TeI < TDMI < 7"DDI 

which is compatible with the fact that the correspond- 
ing couplings are ordered in the following way 



Aei > A 



DMI 



> A 



DDI 



and that r 1/A. 

In a pure material, i.e. without too many impuri- 
ties, it turns out that the fastest recovering of the 
magnetic state and thereby that of the system mag- 
netization occurs via the exchange coupling. In this 
work, we provide details of how this switching occurs. 

We have already started a few experiments for in- 
vestigating the dynamics of coupled thin films grown 



by our collaborators. We intend to perform various 
measurements of FMR with varying field and fre- 
quency using a network analyzer. In addition, the 
slow dynamics of the dimer may be probed by mea- 
suring the isothermal and thermoremanent magne- 
tization by a commercial SQUID in a wide range of 
temperature. 

For a closer comparison with experiments we need 
to consider more general situations with arbitrary di- 
rections of the two anisotropy axes in an oblique mag- 
netic field. Such calculations will be performed nu- 
merically in a subsequent work. 



Appendix A: Details for WC of DDI-MD with MA 

The stationary states are given by 

±a arccos x+ ±b arccos X- 
61,62 = ^ 

±d arccos x+ ±e arccos X- 
0i,02 = ±c7r H (Al) 



where 



xX = ±1, xt = ±1, 



^ see + 9^^ - ^ 9^ + 4^^ 
Xj_ = ±\l — —,x_=± 



36^2 + 16 



16 + 4C 



2 ■ 



and the subindex a...e indicates independence be- 
tween different ±, =f signs. 

In the weak coupling, the minima of the MD are 
located at 



arccos x,— arccos x 
Oi = ± ^ = 



92 = ± 



arccos x_^_+ arccos x_ 



and 



arccos x,— arccos x 
Oi = TTTT ^ = 



±7r =F 



arccos x_^+ arccos x_ 



where the sign of a;_|_ is taken as positive. The maxima 
are at 

arccos x7 — arccos x~ 
61 = ± ^ 



arccos a;, -|- arccos x_ 



62 = ± ± 



and 



arccos x, + arccos x 

61 = iTTT ^ = 

arccos a; 7 — arccos x~ 
O2 = ±TTT = 
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where the sign of is taken as negative. Hence, the 
saddle point are at 

(^1,^2) = (0, ±7r) , (±7r, 0) , (±^, , 



The switching rates for the first and second step in 
the WC regime read 



where 



-(1) 

MAWC 



.(2) 
MAWC 



27r V + 2r) \^ - 2r\{l + ^ + r)\ {I + ^ - r) 



„a+2a£+£^(0) 



47r V (2$ + rp) (1 + 2$ + r) U |2^ - rp\ (1 + 2^ - r^) 



(A2) 



(Cpi) - ipl-i) 



(C^i) -(^^^2) ) +\{W, + 2Qi) , 



with 



i = 1,2 refers to the i* layer and j = 1, 2 refers to the 
j*^ minimum. 
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